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Discovery of kilogauss magnetic fields in three DA white dwarfs 
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Abstract. We have detected longitudinal magnetic fields between 2 and 4kG in three (WD 0446-790, WD 1105-048, 
WD 2359-434) out of a sample of 12 normal DA white dwarfs by using optical spectropolarimetry done with the VLT Antu 
8 m telescope equipped with FORS1. With the exception of 40 Eri B (4kG) these are the first positive detections of magnetic 
fields in white dwarfs below 30 kG. Although suspected, it was not clear whether a significant fraction of white dwarfs contain 
magnetic fields at this level. These fields may be explained as fossil relics from magnetic fields in the main-sequence progeni- 
tors considerably enhanced by magnetic flux conservation during the shrinkage of the core. A detection rate of 25 % (3/12) may 
indicate now for the first time that a substantial fraction of white dwarfs have a weak magnetic field. This result, if confirmed 
by future observations, would form a cornerstone for our understanding on the evolution of stellar magnetic fields. 
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1. Introduction 

The major goal of this work is a better understanding of the 
role played by magnetic fields in the formation and evolution 
of stars. Magnetic fields are already an important ingredient 
during the collapse and fragmentation of protostellar clouds 
which ultimately determines the initial field of pre-main se- 
quence stars. On the main sequence and at later stages of evo- 
lution the magnetic fields have a major impact on the angu- 
lar momentum loss and stellar winds, on building-up chemi- 
cal anomalies and abundance inhomogeneities across the stel- 
lar surface, on convection and the related coronal activity, and 
other evolutionary processes especially in interacting binaries. 
Very strong magnetic fields are d etected in sever al white dwarfs 
and always present in pulsars jMestelll200TI) . The first de- 
tecti on of a magn etic field on a white dwarf was made by 
lKempetal.Nl970l) on Grw+70° 8247, and large spectroscopic 
and polarimetric s urveys have been carried out in the last 
two decades (e.g. iHagen et al. Ill987t iReimers et al 111 9941 
ISchmidt & Smithl M 99.4 IPntnev I M 99.4 iKawka et al. I l2003t> . 
For the white dwarfs the magnetic fields could simply be "fos- 
sil" remnants of the fields already present in main-sequence 
stars, but strongly amplified by contraction. This hypothesis as- 
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sumes that the magnetic flux (e.g. through the magnetic equa- 
tor) is conserved to a large extent during the stellar evolution. 

According to the field amplification theory, the white 
dwarfs play an important role in the investigation of stellar 
magnetic fields. In main-sequence stars magnetic fields have 
been detected directly mainly for peculiar magnetic Ap and Bp 
stars with rather well organized fields and field strengths of the 
order 10 2 - 10 4 Gauss. For weak fields in A to O stars (B < 10 2 
G) direct magnetic field detections are still v ery rare (e.g. 
the field detections in early B stars reported by iNeiner et al. I 
l2003albl) . For sun-like stars ample evidence (coronal activ- 
ity) for the presence of complicated small-scale fields exists, 
but di rect measurements are only possible for the more active 
stars dSaar I [l 9961 iRiiedi et al. Ill997t IValenti & Johns-Krulll 
I2001I) . The contraction to a white dwarf amplifies the magnetic 
fields by about 4 orders of magnitude, so that weak and of- 
ten undetectable magnetic fields on the main sequence become 
measurable during the white dwarf phase. This is supported 
by the known magnetic white dwarfs with megagauss fields 
(B = 10 6 - 10 9 G). Their frequency and space distribution, 
as well as their mass, are consistent with the widely accepted 
view that they are t he descendents of the magnetic Ap and Bp 
stars (e.g. Mathvs 120011) . Another origin seems to be required 
for the magnetic degenerates with weaker fields (unless mag- 
netic flux is lost during the contraction phase). Magnetic main- 
sequence stars with weaker magnetic fields h ave been sug- 
gested as their possible p rogenitor candidates dSchmidt et al. I 
2003t lKawka et al. II2003I) . The B stars on which weaker fields 
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have been detected may be the missing stars. However, even the 
most sensitive observations are limited to some tens of gauss on 
main-sequence stars. 

Thus, magnetic field amplification during stellar evolu- 
tion may offer the opportunity to investigate ~ 1 G magnetic 
fields (averaged global fields) in normal main-sequence stars 
with observations of ~ 1 kG magnetic fields during the white 
dwarf stage. White dwarfs with magnetic fields below lOOkG 
have been either found b y searching for circular polarization 
JSchmidt & Smithll 19941) or by looking for Zeeman splitting 
of n arrow NLTE line co res in the Balmer lines, particularly in 
Ha jKoester et al. Ill998i) . However, the splitting becomes un- 
detectable in intensity spectra for weak fields (< 20 kG) or for 
objects without narrow line core. Therefore, spectropolarime- 
try is the most promising technique for successful detections of 
weak magnetic fields. Up to now detections of magnetic fields 
below 30 kG have not been achieved, exc ept for the very bright 
white dwarf 40 Eri B (V = 8.5), in which lFabrika et al. I §003) 
have detected a magnetic field as low as 4kG. The magnetic 
field detection limit can now be pushed down to a few kG for 
many white dwarfs with spectropolarimetry using 8-10 m class 
telescopes. 

In this paper we present and analyse VLT spectropolari- 
metric data of a sample of 12 white dwarfs in a search for weak 
magnetic fields. In Sect. 2 the observations and data reduction 
are described, while in Sect. 3 the observational method for ob- 
taining the Stokes parameter (V/I) is described. Sect. 4 presents 
the method for determining weak magnetic fields analysing the 
circular polarisation due to a given magnetic field. In that sec- 
tion we also present the results of our analysis, along with the 
description of the x 1 -minimisation procedure applied to our 
data. The determination of the atmospheric and stellar param- 
eters is presented in Sect. 5 and compared with those found 
in the literature. A discussion and conclusions are presented in 
Sect. 6. 

2. Observations and data reduction 

Spectropolarimetric observations of a sample of 12 normal DA 
white dwarfs were carried out during the period 4 November 
2002 - 3 March 2003, in service mode by ESO staff mem- 
bers using FORS1 at the 8 m Unit Telescope 1 (UT1) of the 
Very Large Telescope, 'Antu". FORS1 is a multi-mode fo- 
cal reducer imager and grism spectrograph equipped with a 
Wollaston prism and rotatable retarder plate mosaics in the par- 
allel beam allo wing linear and ci rcular polarimetry and spec- 
tropolarimetry ( App enzeller et alIT9 98). 

With a 0.8" wide slit we obtained a (FWHM) spectral res- 
olution of 4.5 A. The data were recorded with a backside- 
illuminated thinned, AR coated Tek. 2048x2048 CCD with 24 
yum pixels which correspond to a pixel scale of 0.2"/pixel in 
spatial and 1 A/pixel in spectral direction. 

Spectra were acquired with grism G600B, which allows ob- 
servations in the spectral range 3400-6000 A, covering all H i 
Balmer lines from Hfi to the Balmer jump simultaneously. A 
reflex image from the FORS optics affects the wavelength re- 
gion from 4000 to 4100 A which corresponds to the blue wing 
of the H6 line. Although the reflex shows up in the intensity 



spectrum it seems to produce no spurious signal in circular po- 
larisation. The reflex is known to occur for this pa rticular grism 
G600B / Wollaston configuration (e.g. ISchmid et ai~l2~0 03). 

The observations were split into several cycles to avoid 
saturation. The feasibility of circular spectropolarimetry with 
the required high signal -to-noise ratio u sing FORS 1 had been 
demonstrated by Bagnulo et al. ] d2002l |2004) with test mea- 
surements of one magnetic and one non-magnetic A-star. 
Stellar rotation may cause continuous changes in the field 
orientation and, therefore, in the polarisation signal. Hence, 
we have split the observations of an individual target into 
more than one polarisation measurement taken during different 
nights (see Sect. 3). In this way, there is a much lower probabil- 
ity that a candidate with a magnetic field escapes detection due 
to a special orientation at the time of the observations (where 
circular polarisation cancels). Hence, having for most objects 
data for more than one epoch, we are able to assess for the 
presence of rotational modulation of a possible magnetic field. 

2.1. Selection of the sample 

We have searched for the circular polarisation signatures of the 
Zeeman effect in the Balmer lines in 12 of the brightest south- 
ern white dwarfs (having 11.4 mag < V < 14 mag) with RA 
between h and 12 h . 

Our sample was carefully selected on the basis of VLT- 
UVE S spectra taken within the SPY survey jNapiwotzki et al. I 
l2003h : SPY (Supernovae type la Progenitor surveY) is a ra- 
dial velocity search for close binary systems of two white 
dwarfs (double degenerates; DD). If these systems are close 
enough they will merge due to gravitational wave radia- 
tion and if the combined mass of these mergers exceed the 
Chandrasekhar limit for white dwarfs these are potential pro- 
genitors of Supernovae type la. 

SPY was carried out with the high-resolution Echelle spec- 
trograph UVES at the Kueyen (UT2) of VLT. With the set-up 
used for SPY UVES provides almost complete spectral cov- 
erage of the wavelength range from 3200 A to 6650 A, with 
a spectral resolution of 0.3 A. For more details, please refer to 
Nap iwotzki et al. IfcOOliEoOsT) . SPY observations were used to 
check candidates for our project for spectral peculiarities and 
magnetic fields strong enough to be detected in intensity spec- 
tra. Hence, our targets were selected with the criterion of not 
having any sign of Zeeman splitting visible in the SPY spectra, 
and hence no magnetic fields above a level of about 20 kG. 

All our targets have strong hydrogen lines, ideal for mea- 
suring line polarisation, and no peculiarities (such as MG- 
magnetic fields or a bright companion or any indication of mag- 
netic fields). 

2.2. Data reduction 

Data reduction was performed using midas routines. Flat field, 
bias and wavelength calibration exposures were taken during 
day time with the same instrument setup, after each observing 
night. 
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In all frames the bias level was subtracted and the frames 
were cleaned of cosmic ray hits. For each observing night a 
unique master flat-field was calculated from the median flux 
of all flat-fields taken that night. The stellar spectra were ex- 
tracted from the flat-field corrected frames as a sum over about 
forty CCD rows for each beam. Sky spectra were obtained from 
adjacent regions (about twenty CCD rows) on the detector to 
the observed stars (below and above) and subtracted from the 
object spectra. Wavelength calibration was done using HgCd, 
He and Ar arc spectra, which was independently performed for 
wavelengths of the ordinary and extra-ordinary beams with the 
corresponding beams of a reference spectrum. 

Stokes / was obtained as a sum over all beams, while the 
calculation of Stokes (VI I) is described in the following sec- 
tion. 

3. Circular Polarisation 

In order to obtain circular polarisation spectra a Wollaston 
prism and a quarter-wave plate were inserted into the optical 
path. Each exposure yields two spectra on the detector, one 
from the extra-ordinary beam and the other from the ordinary 
beam. 

In the ideal case, the polarisation information (one Stokes 
parameter per exposure) is contained in the ratio, at each wave- 
length, of the intensities in the two spectra (from the ordinary 
and extra-ordinary beams) but it is mixed up with the system 
gain ratio for the pixels concerned. Alternatively, the effect of 
the not well defined pixel gain can be significantly reduced by 
inverting the sign of the polarisation effects in a second expo- 
sure (by rotatin g the /l/4-plate by 90°), whi le leaving the gain 
ratios identical (Tinber gen & Ruttenll 19971) . Hence, Stokes V 
is obtained from a differential measurement of photon counts 
in either the ordinary or extra-ordinary beams, measured at two 
different angles of the retarder waveplate. In this way errors 
from changes in the sky transparency, atmospheric scintilla- 
tion, and various instrumental effects are significantly reduced, 
so that photon noise remains as the dominant error source. 

Another important source of sytematic error when obtain- 
ing Stokes V could be the wavelength calibration procedure. 
An incorrect calibrated wavelength scale for each of the two 
analysed spectra will yield the line profiles of the two spectra 
not to be perfectly aligned, even in the absence of magnetic 
field, leading to spurious polarization signals in every line, that 
change sign as the wave plate rotates. No such spurious signals 
were detected during the calibration process, confirming that 
FORS is a very stable instrument. Even if spurious line signals 
at the noise level are present, they would be compensated at 
least partly by the combination of data taken with retarder plate 
position angles +45° and -45°. We note that wavelength cali- 
brations were made for each observing date separately. Because 
most targets were observed for two different dates we could 
check that no spurious magnetic field detection due to wave- 
length calibration errors are present. 

We adopted the FORS1 standard observing sequence for 
circular polarimetry consisting of exposures with retarder plate 
position angles +45° and -45°. The number of exposures, n, 
is reported in Table 1 for each object. For instance, for n = 4 



the sequence of position angles would go: +45°, +45°, -45°, 
-45°. 

In order to derive the circular polarisation from a sequence 
of exposures, we added up the exposures with the same quarter- 
wave plate position angle. The Stokes (V/I) can be obtained as 



V = (R-l) 

i (R + iy 



withtf 2 = i4 



\Jo/ a= _ 45 



(1) 



where V is the Stokes parameter which describes the net cir- 
cular polarisation, / is the unpolarized intensity, a indicates 
the nominal value of the position angle of the retarder wave- 
plate, and f and f e are the fluxes on the detector from the 
ordinary and extra-ordinary beams of the Wollaston, respec- 
tively. Eq.[T]is equi valent to formula 4 .1 in the FORS 1 + 2 User 
Manual (Szeifert & Bohnhardt"1l2003l) . In lDonati et all i 19971) 
a detailed description of the method for the suppression of spu- 
rious signals in spectropolarimetric observations that we ap- 
plied is provided. 

The 12 selected white dwarfs are listed in Table 1 together 
with some characteristics of the observations. The Heliocentric 
Julian Dates correspond to the beginning of each observing se- 
quence. 

4. Determination of weak magnetic fields 

4.1. Theoretical polarisation 

For the weak magnetic fields expected in our white dwarfs 
the longitudinal Zeeman effect (due to magnetic fields paral- 
lel to the line of sight) is measurable for the broad Balmer 
lines with intermediate spectral resolution. This is possible be- 
cause, in first approximation, the positive and negative extrema 
of Stokes V occur at the wavelengths where the second deriva- 
tive of Stokes I changes its sign. These points could possibly 
be regarded as the transitions between the core and the wings 
of the line. Hence, the wavelength difference between the max- 
ima of the opposite net circular polarisation in the two wings of 
the line, is correlated with the width of the core of the line and 
not with the Zeeman splitting AAr , which can be much smaller. 
The detection threshold therefore depends mainly on the noise 
level in V and the width of the core of the spectral lines, and far 
less on the achieved spectral resolution. 

With FORS 1 on the VLT we reach a noise level in circular 
spectropolarimetry of about 3 x 1(T 3 / C , where I c is the con- 
tinuum intensity. This allows us to detect magnetic fields of a 
few kG for the brightest white dwarfs in our sample with ex- 
posure times of about 1 hour. T his is a s ignificant improvement 
compared to the observations of lSchmidt & Smith I (ll995). who 
reached a detection limit of about 20 kG (2.3<x limit) with a 2 m 
class telescope. By adding the signal of several Balmer lines, 
which are observed simultaneously with our instrument setup, 
the magnetic sensitivity can be further enhanced. However, 
the error of the magnetic field determination increases for the 
higher series numbers, so that a reliable analysis can be based 
on H/3 and Hy only. 

The theory of spectral line formation in a magnetic atmo- 
sphere shows that the splitting of a spectral line observed in 
both senses of circular polarisation is proportional to (B z ), the 
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Table 1. Details of VLT observations. T he provided a and 6 coordinates refer to epoch 2000 as measured in the course of 
the SPY project (see iKoester et al, 1 12001 ). Spectral types, r sp , and measured V magnitudes were taken from the catalogue of 
lMcCook&Sion1 l ll999l) . 



Target 


Alias 


a 


5 


V 


HJD 


*exp 


;•; 












(mag) 


(+2452500) 


(s) 






WD 0135-052 


LHS 1270 


01 37 59.4 


-04 59 45 


12.84 


108.584 


546 


4 


DA7 


WD 0227+050 


Feige 22 


02 30 16.6 


+05 15 51 


12.65 


137.601 


381 


6 


DA3 












169.545 


381 


6 




WD 0310-688 


CPD-69 177 


03 10 31. 


-68 36 03 


11.40 


195.536 


119 


14 


DA3 


WD 0346-011 


GD50 


03 48 50.2 


-00 58 31 


13.99 


137.657 


699 


4 


DAI 












174.600 


699 


4 




WD 0446-789 


BPM 03523 


04 43 46.4 


-78 51 50 


13.47 


109.710 


699 


4 


DA3 












168.569 


699 


4 




WD 0612+177 


T rTTTI 1 1 O 1 O 

LTT 1 1818 


06 15 19.0 


+ 17 43 48 


13.39 


109.756 


699 


4 


DA2 












172.561 


699 


4 




WD 0631 + 107 


KPD063 1 + 1043 


06 33 50.6 


+ 10 41 09 


13.82 


200.612 


699 


4 


DA2 












202.607 


699 


4 




WD 0839-327 


CD-32 5613 


08 41 32.6 


-32 56 35 


11.90 


108.803 


186 


10 


DA6 


W U Uo Jy — U Jy 




HQ m 17 *\ 
Uy UZ 1 / . J 


CIA H7 1 7 
— U^f U / 1Z 


1 9 An 


1 1A 7HQ 


Z07 


Q 
o 


1 A \ _ 












196.701 


289 


8 




WD 1042-690 


BPM 06502 


10 44 10.5 


-69 18 20 


13.09 


168.832 


426 


6 


DA3 












174.755 


426 


6 














195.783 


426 


6 




WD 1105-048 


LP 672-001 


1 1 08 00.0 


-05 09 26 


12.92 


141.838 


486 


4 


DA3 












169.791 


486 


4 




WD 2359-434 


L 362-81 


00 02 10.7 


-43 09 55 


13.05 


83.510 


546 


4 


DA5 












108.542 


546 


4 





average of the component of the magnetic field along the line 
of sight averaged over the visibl e stellar he misphere, i.e., the 
mean longitudinal magnetic field llBabcockll 947. 1958). 

For field strengths below lOkG the Zeeman splitting of the 
Balmer lines is less than approximately 0.1 A. This is well be- 
low the width of the cores of the Balmer lines in all the stars 
of our sample (typically a few A ). Therefore, we can apply 
the weak-field approximation (e.g..lAngel & Landstreetll970t 
lLandi deglTnnocenti & Land i degl'Innoce ntil Il973h without 
any loss of accuracy. According to this approximation the mea- 
sured V and / profiles are related to (B z ) by the expression: 

7 - -SeffQA 2 --^) , (2) 

/ I OA 

where g e fj is the effective Lande factor (= 1 for all hyd rogen 
lines of any series. Icasini & Landi deglTnnocenti1 [l994). A is 
the wavelength expressed in A, (B z ) is the mean longitudi- 
nal component of the magnetic field expressed in Gauss and 
the constant C z = e/(4mn e c 2 ) (^ 4.67 x 10~ 13 G~' A~'). Note 
that this approximation also holds if instrumental broadening 
is present, but it is not generally corre ct if the profiles are rota- 
tionally broadened ( Landstree t 1 1 982h . 

Eq.|3]is evaluated by using the high signal-to-noise spectra 
I{A) and their derivative dl/dA. For a given wavelength A, we 
approximate dl/dA by the average of - 7,)/(/l, + i - A,) and 
(Ii-Ii-^/tti-Ai-i). 

The error associated with the determination of the longitu- 
dinal field obtained from individual Balmer lines is larger for 
Balmer lines at shorter wavelengths than for lines at longer 



wavelengths. This is due to the combination of two effects: 
while the Zeeman effect increases as lambda squared, most 
other line broadening effects depend linearly on lambda, so that 
the magnetic field is better detected at longer wavelengths than 
at shorter wavelengths; furthermore, the Balmer lines at shorter 
wavelengths are less deep, so that dl/dA is smaller. Using H/3 
and Hy simultaneously, we obtained a determination of the 
mean longitudinal magnetic field that best fit the observed (V/I) 
(as explained in the next section). 

4.2. Fitting of the magnetic field 

In Fig.^we present the normalised spectra of our sample of 12 
DA white dwarfs. Multiple observations of the same star were 
added up in order to achieve a higher signal-to-noise level. In 
some spectra, particularly prominent in WD 0135-052, a small 
bump is visible on the blue side of the H6 line, due to the reflex 
mentioned in Sect. 2. 

For the circular polarisation spectra {V/I) flux weighted 
means were calculated according to the following formula for 
two measurements 1 and 2: 

\iL h+h ' ' 

For the averaging we have checked that the effect of the dif- 
ferent barycentric corrections for multiple observations of the 
same object is much smaller than our spectral resolution. The 
averaging helps to extract signal hidden in the noise of the 
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3800 4000 4200 4400 4600 

Wavelength (A) 



4800 



5000 



Fig. 1. Normalised spectra of our sample of white dwarfs in the region of the Balmer series. When more than one sequence of 
exposures was available for a white dwarf, the average spectrum is presented. The positions of the Balmer lines are indicated. 
The spectra are displaced vertically by 0.7 units (or multiples thereof) for a better visualisation. 



individual exposures if the stellar Zeeman signal remains un- 
changed with time. 

Fig. |2 shows the circular polarisation spectra (V/I) of our 
sample of white dwarfs in the region of the Balmer series (no 
binning along the spectrum was applied to the data). The po- 
sitions of the Balmer lines are indicated by vertical dashed 
lines. In two of the averaged spectra (corresponding to stars 
WD 0446-789 and WD 2359-434) a moderate S-shape circu- 
lar polarisation signature across HB and Hy can be noticed. For 
WD 1105-048 polarisation reversals of those lines were only 
present in the observation of 9 Jan. 2003 (see Fig.|5|l. 

In order to determine the longitudinal component of the 
magnetic field for each measurement we compared the ob- 
served circular polarisation for an interval of +20 A around HB 



and Hy with the prediction of Eq.|2] The best fit for (B : ), the 
only free parameter, was found by a ^-minimisation proce- 
dure. If we assume that no magnetic field is present, all de- 
viations from zero polarization are due to noise. This can be 
expressed by the standard deviation <x over the respective in- 
tervals around the Balmer lines. If the best fit magnetic field 
is indeed close to zero, the reduced y 2 should be automati- 
cally close to unity. However, in the case of a finite magnetic 
field not all deviations from zero would be due to noise, so that 
the best fit would h ave a reduced x 1 smaller than 1 . Following 
Pres s et al. I ll 19861) we determined the statistical error from the 
rms deviation of the observed circular polarisation from the 
best-fit model. The lcr (68.3%) confidence range for a degree 
of freedom of 1 is the interval of B z where the deviation from 
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I | i i i I I i i wd2359-434 « , i 



0.10 



0.08 




C 0.06 
> 



0.04 



0.02 



0.00 



wdl 105-048 



I I i I i i 

i i wdl()42-690 



wd0359-039 



I I I I I I I I I . i 



wd0839-327 



i 

wd06 




wd0631+107 



wd0612+177 



wd0446-789 




wd0346-01,l 



wd0310-688 



wd0227+050 



wd0135-052 



ii. n 



i i i j i r ^ m> I 



3800 4000 4200 4400 4600 

Wavelength (A) 



4800 



5000 



Fig. 2. Circular polarisation spectra (V/I) of our sample of white dwarfs in the region of the Balmer series. The average (V/I)- 
spectrum is plotted for multiple observations. For all spectra a horizontal line is drawn to indicate the zero level. Dashed vertical 
lines represent the positions where the Balmer lines are located. Spectra are displaced with relative shifts for a better visualisation. 



the minimum is Ajf 2 = 1. Note that this is a purely statistical 
error and does not account for systematic errors, e.g. from the 
limitation of our low-field approximation which, however, we 
expect to be rather small, or from wavelength calibration errors. 

The results for each of our single observations is sum- 
marized in Table[2] where we have listed the best fits for B, 
for the HB and Hy lines as well as the weighted value B z = 
(B^yWy + B Z fiWp)l(w y + w p) where Wi = 1/of (i = y,B). The 
probable error is given by <x = (2 WiY 1 ^ 2 . When multiple ob- 
servations were added up to reduce the noise level (by means 
of Eq.0, results are labeled as average. 

In three stars we find a significant magnetic field: 
WD 0446-789 (4283 + 640 G, see Fig©, WD 2359-434 
(-3138 + 422G, see Fig.EJ, and WD 1105-048 (-2134 + 



447 G). Positive and negative signs of B z indicate opposite 
magnetic polarity. For the first two stars the magnetic field is 
detected at the 3<r level individually from HB and Hy (at least 
for the averaged spectra), as well as from the combination of 
both lines. This increases our confidence in both detections (cf. 
Sect. 4.3). 

Although several lines of the Balmer series are available, 
the higher members of the series do not contain enough V- 
signal to give reliable results. However, the analysis of H6 (the 
only line among the higher members with its V-signal exceed- 
ing the l<x(V)), confirms the positive detections of magnetic 
fields in the three objects WD 0446-789, WD 1105-048 and 
WD 2359-434. 
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Table 2. Magnetic fields derived from the Hy and HB lines for our sample of white dwarfs. B(cr) provides the magnetic field in 
units of the cr level. Detections exceeding the 3cr levels are given in bold. cr(V) is the standard deviation of the observed {VI I)- 
spectrum obtained in the region 450CM-700 A. Lower limits on the detectability of the magnetic field from the line polarisation 
peaks, calculated at the lcr level of the noise, are given in the last two columns. Multiple observations that were averaged prior 
to analysis are labeled average. 



Target 


Date 


<r(V) 


B(G) 




S(G) 


fl(o-) 


B(G) at lcr 






(io- 3 / c ) 


Hy 




Hy,/5 




Hy 




WD 0135-052 


29/11/02 


1.2 


1112 ±625 


492 ± 373 


654 ± 320 


2.04 


1277 


861 


WD 0227+050 


28/12/02 


1.0 


493 ± 1177 


-776 ± 707 


-439 ± 606 


0.72 


2486 


-1928 




29/01/03 


1.0 


-710 ±971 


1854 ±709 


962 ± 572 


1.68 


-2969 


1989 




AVERAGE 


0.7 


-105 ±789 


590 ± 469 


408 ± 403 


1.01 


-1924 


1194 


WD 0310-688 


24/02/03 


0.8 


-504 ± 741 


114 ±548 


-104 ±440 


0.24 


-2106 


1404 


WD 0346-011 


28/12/02 


1.8 


-9299 ±8510 - 


-6815 ±5295 


-7508 ± 4495 


1.67 


-36860 


-16600 




03/02/03 


2.6 


3222 ± 9374 


7542 ±6717 


6076 ± 5459 


1.11 


28440 


24820 




AVERAGE 


1.1 


-5967 ± 7087 - 


•1496 ±4226 


-2668 ± 3629 


0.74 


-23060 


-10680 


WD 0446-789 


30/11/02 


1.3 


417 ± 1403 


4251 ±910 


3115±763 


4.08 


3731 


2340 




28/01/03 


1.4 


7771 ± 1476 


5370 ± 1196 


6321±929 


6.80 


4118 


2815 




AVERAGE 


0.9 


3612 ± 1005 


4743 ± 832 


4283±640 


6.69 


2495 


1594 


WD 0612+177 


30/11/02 


1.3 


2416 ± 1456 - 


•1148 ± 1203 


297 ± 927 


0.32 


4060 


-2620 




01/02/03 


1.4 


86 ± 1154 


236 ± 1171 


156 ±821 


0.19 


3742 


2915 




AVERAGE 


0.9 


1132 ±920 


-457 ± 789 


216 ±598 


0.36 


2672 


-1939 


WD 0631 + 107 


01/03/03 


2.2 


-1310 ± 1785 


-170 ±1815 


-749 ± 1272 


0.59 


-7300 


-5890 




03/03/03 


2.1 


-790 ± 1870 


310 ± 1775 


-211 ± 1287 


0.16 


-6630 


4340 




AVERAGE 


1.2 


-990 ± 1440 


260 ± 1090 


-195 ±869 


0.22 


-4010 


2630 


WD 0839-327 


29/11/02 


0.7 


-313 ±325 


402 ± 261 


121 ± 203 


0.60 


-1008 


761 


WD 0859-039 


03/02/03 


1.5 


2083 ±1164 


11 ±987 


877 ± 752 


1.17 


4400 


2727 




25/02/03 


1.2 


2183 ± 1372 


-401 ±911 


389 ±758 


0.51 


3234 


-2222 




AVERAGE 


0.8 


1186 ±928 


382 ± 607 


622 ± 507 


1.23 


2319 


1511 


WD 1042-690 


28/01/03 


1.3 


-719 ±2333 - 


•2850 ± 2301 


-1799 ± 1638 


1.10 


-3624 


-3247 




03/02/03 


1.0 


1627 ± 1464 


-865 ± 1163 


99 ±910 


0.11 


2963 


-3560 




24/02/03 


0.9 


-2374 ± 1060 


103 ± 864 


-885 ± 669 


1.32 


-2764 


1804 




AVERAGE 


0.6 


141 ± 745 


-406 ± 620 


-182 ±476 


0.38 


1833 


-1508 


WD 1105-048 


01/01/03 


1.1 


-1153 ±632 


55 ± 807 


-693 ± 497 


1.39 


-2863 


1599 




29/01/03 


1.1 


-4686 ± 1032 


-3305 ± 979 


-3959±710 


5.58 


-2756 


-2164 




AVERAGE 


0.8 


-2760 ± 598 


-1337 ±675 


-2134±447 


4.77 


-2156 


-1538 


WD 2359-434 


04/11/02 


1.0 


-3003 ± 1433 - 


-5719 ± 1289 


-4504±958 


4.70 


-2914 


-2975 




29/11/02 


1.0 


-1905 ±702 


-4224 ± 796 


-2919±526 


5.55 


-4880 


-3907 




AVERAGE 


1.0 


-2066 ± 546 


-4738 ± 667 


-3138±422 


7.44 


-1785 


-1532 



In WD 0446-789 the analysis of the 30 Nov. 2002 observa- 
tion around Hy yields a magnetic field of 417 + 1403 G, corre- 
sponding to a xtm n = 1-0 (instead of x 2 min - 0.6 for the 28 Jan. 
2003 observation), i.e. a O.lcr level detection. If we do not use 
this outlier for our combined result from the two observations, 
we obtain a 3cr detection of 4743 + 832 G. 

In the case of WD 1105-048 the different values of the 
magnetic field obtained for the two observations on 1 Jan. 2003 
(-693 ± 497 G) and 29 Jan. 2003 (-3959 + 710G, see FigJJi 
differ significantly. This may indicate a different orientation of 
the magnetic field between the two epochs due to stellar rota- 
tion. 

The case of WD 0135-052, where the formal value for the 
combined result from HB and Hy exceeds the lcr error by a 
factor of two, is rather uncertain since it is based on one ob- 
servation only. In addition, this object is a double-lined binary 
with a period of 1 .556 days and it is not clear whether the or- 



bital variation can mimic a magnetic field. We agree with the 
referee, who pointed out that this could just be a 2cr level detec- 
tion, which can be expected for our sample size of field mea- 
surements just by chance. However, WD 0135-052 is a good 
candidate for follow-up observations, although one must bear 
in mind that the magnetic field would be diluted by the pres- 
ence of a non-magnetic companion. 

In all other stars in our sample the best-fit magnetic field is 
below or close to the lcr error. 

In Figs. [5] to |6] the observed circular polarisation spectra 
of the white dwarfs (thin solid line) are presented together 
with the {VI I) spectra predicted by the low -field approxima- 
tion (thick solid line, red colour in the electronic version) for 
the given value of B. The light shading (cyan in the electronic 
version) shows the values of (V/I) obtained varying B by the 
indicated errors, and can be more clearly distinguished from 
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B = 4283±640 G 
(V/I) 




B = -3138±422 G 
(V/D 





-0.004 
-0.006 

4200 4250 



0.006 r 

0.004 r 

0.002 

0.000 
-0.002 
-0.004 
-0.006 : 
4200 



B = -313R±422 G 
(V/I) 











-0 

-0 
-0 




0.006 F 



B = -3138±422 G 
(V/I) 



4000 
Wavelength (A) 



Fig. 3. Circular polarisation spectra (V/I) of WD 0446-789 
(average of observations from 30/1 1/03 and 28/01/03, thin solid 
line) in the region of H/3 (top), Hy (middle) and close to the 
Balmer jump (bottom panel). The solid horizontal line indi- 
cates the zero level. The horizonal dashed lines indicate the po- 
sition of the lcr level of the average (V/I) spectrum. The light 
shading (cyan in the electronic version) represents the variation 
between (V/I) spectra predicted by the low-field approximation 
(Eq.|2j using magnetic field values of 4283+640 G (thick solid 
line, red in the electronic version). 



the predicted (V/I) spectra at the peaks of the S-wave features 
(in particular in the region of Hy). 

4.3. Accuracy of the fitting procedure 

Although the x 1 fitting provides formal errors, it is difficult to 
judge whether this error actually provides a reasonable esti- 
mate of the confidence range for our purpose. For example, 
the choice of our interval of ±20 A around H/3 and Hy is arbi- 
trary. Therefore, we performed simulations for different mag- 
netic fields, interval sizes, and noise levels in order to test our 
approach. The Gaussian noise for our simulated polarisation 
spectra was provided by the random number generator imple- 
mented in IDL. 

As an example we used the solution for the average of both 
WD 0446-789 observations (B=4283+640 G), its noise level 
(cr = 0.001 at both H/3 and Hy) and calculated 1000 simulated 




3900 4000 

Wavelength (A) 



4200 



Fig. 4. As in Fig.|3]but for WD 2359-434 (average of obser- 
vations from 04/11/03 and 29/11/03), where the best fit is at 
(B z ) = -3138 + 422 G. 



polarisation spectra for the interval size of +20 A. The aver- 
age result for the magnetic field was: 4301+442 G (solution 
for Hy,/3), 4276+554 G (for H/3 only), and 4339+725 G (for 
Hy only). The mean value of B agrees very well with the pre- 
scribed value. The mean of the individual standard deviations 
of the 1 000 runs are 4 1 3 G, 54 1 G and 667 G, for Hy , /3, H/3 and 
Hy, respectively. This is slightly smaller than 640 G, 832 G and 
1005 G from the WD 0446-789 observation. This may be due 
to the fact that we assume a homogeneous magnetic field which 
may not be the case in reality, or it may indicate a small hidden 
source of errors in the data (e.g. a non-Gaussian distribution of 
the noise). 

Fig-Elshows one of the artificial polarisation measurements 
from our 1000 simulations. The fit is 4449+449 G in this partic- 
ular example. It demonstrates that a highly significant detection 
is possible at our noise level. 

The resulting magnetic field changed by less than 2% if an 
interval of + 40 A around H/3 and Hy was used. In the case 
of a smaller interval, however, larger deviations occur since not 
the whole range where the predicted (V/I) differs from zero is 
included. 

We are aware of the fact that the polarisation feature is 
mostly buried in noise in the (V//)-spectrum, with the excep- 
tion of the narrow peaks on both sides of the hydrogen line cen- 
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B = -3959±710 G 
(V/I) 




0.006 
0.004 
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-0.004 
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B = 4449 ±449 G 
(V/I) 

J, 
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Fig. 5. As in Fig.|3]but for the 09/01/03 observation of the white 
dwarf WD 1105-048 where the best fit is at (B z ) = -3959 + 
710G. 
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Fig. 6. Artificial polarisation spectrum with the same noise 
level as WD 0446-789 plotted in Fig.|3] For this particular ex- 
ample, out of the 1000 simulations, the fitting procedure results 
in (B z ) = 4449 + 449 G. 



ters. Therefore, on a plot the magnetic information is, to a large 
extent, invisible to the eye since in the wings the solution is 
based on an average small excess of right- and left-handed po- 
larisation on different side of the line core, respectively, which 
clearly contributes to our x 2 analysis. However, if one wants 
to obtain an idea of what magnetic field is needed so that the 
predicted polarisation peaks exceed the lcr level of the noise, 
we provide lower limits for this detectability in columns 8 and 
9 of Table|2] The cr was obtained as the standard deviation of 
the observed (V//)-spectrum in a region with no lines over 200 
A (i.e. 4500-4700 A), for each single observation. 

We note that in the three objects (WD 0446-789, 
WD 1105-048 and WD 2359-434) with a positive detection 
according to the x 1 analysis, the peaks of the S-wave circu- 
lar polarisation signature reaches the 2<x level of the noise so 
the detection is, therefore, confirmed by this very conserva- 
tive approach. For the averaged spectra of WD 0446-789 and 
WD 2359-434, the detections derived from the more reliable 
HB line even reach the 3<x level. 



5. Atmospheric and stellar parameters for the 
target stars 

For a discussion of the evolutionary status of our white dwarfs 
with and without detected magnetic fields we need to know 
the fundamental parameters of our objects. Masses and cool- 
ing ages are of special interest to distinguish between the pro- 
posed formation scenarios. These quantities can be computed 
from the fundamental stellar parameters temperature and grav- 
ity, which can be derived by a model atmosphere analysis of 
the spectra, and theoretical cooling tracks. 

Since all white dwarfs selected for our project are bright by 
white dwarf standards, at least one model atmosphere analysis 
was published for each object. This is a somewhat inhomoge- 
neous collection of data relying on spectra of different quality 
and a variaty of methods, including analysis of UV spectra and 
parallax measurements. Most of the programme stars have al- 
ready been observed with the SPY project. Although our spec- 
tra have lower resolution than the SPY spectra, which aim at 
the measurement of radial velocity variations seen in the NLTE 
core of Ha, our flux spectra have a very high signal-to-noise 
ratio, allowing a very accurate determination of the effective 
temperatures and gravities. 
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The observed line profiles are fitted with theoretical spectra 
from a large grid of NLTE spectra c alculated with the NLTE 
code developed by IWerner I ( 119861) . Basic assumptions are 
those of static, plane-parallel atmospheres in hydrostatic and 
radiative equilibrium. The adopted chemical composition is 
pure hydrogen, which is appropriate for the stars of our sample. 
A description of th e model ca lculations and the adopted atomic 
physics is given in iNapiwotzki et al. I Jl999l) . where a discus- 
sion of the impact of NLTE on the atmospheres and line pro- 
files of DA white dwarfs is provided as well. The coolest four 
white dwarfs of our sample (WD 0310-688, WD 0839-327, 
WD 1105-048 and WD 2359-434) were analysed with a grid 
of LTE model spectra computed by D. Koester for the analysis 
of DA white dwarfs. The input physics of these model atmo- 
spheres is described in some detail in Finle v et al. 1 dl997t) .LTE 
models are more reliable below 17000 K, because the NLTE 
models do not take into account convection and collision in- 
duced absorption by hydrogen quasi-molecules. On the other 
side, NLTE effects are small at these temperatures (Napiwotzki 
et al. 1999) and of no relevance for our results. 

The line fits were per formed with a mo dified version of the 
least-squares algorithm of Berge ron et al. I d 19921) described in 
INapiwotzki et al. | dl999h . The observed and theoretical B aimer 
line profiles are normalised to a linear continuum in a con- 
sistent manner. Radial velocity offsets are corrected by shift- 
ing the spectra to a common wavelength scale. The synthetic 
spectra are convolved to the resolution of the observed spec- 
tra (4.5 A) with a Gaussian and interpolated to the actual pa- 
rameters. The atmospheric parameters are then determined by 
minimising the x 2 value b y means of a Lev enberg-Marquardt 
steepest descent algorithm dPress et al. Ill986l) . This procedure 
is applied simultaneously to all Balmer lines of one observed 
spectrum. Formal errors can be derived from the covariance 
matrix. However, due to the very high signal-to-noise ratios of 
our spectra, the estimated statistical errors are extremely small 
(not larger than 20 K in temperature and 0.005 dex in the log 
of gravity) and do not provide realistic estimates for the whole 
error budget. The real external errors for the white dwarfs in- 
vestigated in this ar ticle can be estimated t o be 2.3% in T e g and 
0.07 dex in logg (Napiwotzki et al. 1 19991) . 

Our fit results are provided in columns two and three of 
Table|3] The parameters are the straight average of the fit results 
for the individual spectra. For the reasons discussed above, we 
do not give the formal fit errors. WD 0135-052 was not fitted 
by us, because this is a well known d ouble-lined binary con- 
sisting of two DA stars (Safferet al. 1988), which requires a 
special treatment. White dwarf masses were computed from a 
comparison of parameters derived fro m the fit with the grid 
of wh ite dwarf cooling sequences of iBenvenuto & Althaus I 
d 19991) for an envelope hydrogen mass of 10" 4 Mwd, and are 
given in column four of Table [3] A sample fit is shown in 
Fig.0 Most Balmer lines are fitted very well, but the observed 
profile of HS is only poorly reproduced. This results from the 
reflex in the FORS1 optics mentioned in Sect. 2. Thus we ex- 
cluded this line from our fitting procedure. However, parame- 
ter changes are small, if H6 is included anyway. Spectroscopic 
distances c/(spec) are determined from the absolute magnitudes 
computed for the given stellar parameters using the synthetic V 
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Fig. 7. Model atmosphere fit (solid line) of the 30/1 1/02 obser- 
vation of WD 0446-789 (solid histogram). The dashed model 
profile for H6 indicates that this line was not used during the 
fitting procedure. 

band fluxes computed by Berger on et al. I (fl995) and the mea- 
sured V magnitudes given in our Table 1 , and are provided in 
column five of Table Note that the error limits for T e g and 
log g given above correspond to a 7% error in the distance. 

Tablets supplemented by physical parameters of our sam- 
ple of white dwarfs collected from the literature. Fundamental 
parameters such as temperature and gravity determined by 
previous spectroscopic studies are given in columns six and 
seven. The agreement is generally good. Differences are usu- 
ally within the error limits discussed above. Rotational veloci- 
ties v sin i, if available, are shown in column eight. If parallax 
measurements exist, the resulting distances d(trig) are given in 
column nine (H indicates Hipparcos values, while Y are ground 
based data from the Yale general catalogue of trigonometric 
parallaxes). 

5.1. Comments on individual objects 

WD 0135-052 (L 870-2): We have only one single observa- 
tion of this target, which indicates a very weak magnetic field 
of 650 G on a 2<x level (see Sect. 4.2). It is reported as a 
double-lined spectroscopic binary composed of a detached pair 
of DA whit e dwarfs with an o rbital period of the system of 
1.556 days dSafferet al.lll 988). The system h as a mass rati o 
q = Mi IM 2 = 0.90 dSaffer et alll988HBergeron et al. 1(79 89). 
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Table 3. Fitted parameters of the white dwarfs and supplementary data from literature. The three objects with positive detections 
of magnetic fields are labelled in bold face. 



















literature 


WD 


T e ff 




M d(spec) 


^cool 


T c ff 


l°g£ 




v sin; rf(trig) 




(kK) 




/AT \ 


(pc) 


(Myr) 


(kK) 






(kms~') (pc) 


0135-052 






0.47/0.52 1 ' 14 






7.47/6.92 


7.80/7 


eoi.14 


YZ.j + U.4Y 


0227+050 


18.45 


7.79 


0.51 


26.1 


65 


19.07 
18.21 


7.78 2 
7.95 16 




24.3 ± 2.9H 16 


0310-688 


15.73 


8.00 


0.61 


11.2 


170 


15.71 
16.18 


8.16 16 
8.06 4 




0±4 10 10.2 ± 0.2H 16 


0346-011 


43.17 


9.08 


1.24 


34.0 


59 


39.51 
43.20 
40.54 
43.10 


9.07 12 
9.21 18 
9.22 2 
9.09 4 




<5U 


0446-789 


23.45 


7.72 


0.49 


48.8 


21 


23.20 
23.61 


7.70 17 
7.83 4 






0612+177 


25.58 


7.86 


0.56 


46.1 


17 


24.66 
25.17 
25.94 


7.94 2 
7.83 5 
7.97 4 




36.1 ±2.9Y 15 

References: 


0631+107 


27.25 


7.81 


0.54 


61.9 


13 


26.59 
27.20 
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8.00 s 






0839-327 


9.24 


7.90 


0.54 


8.3 


660 


8.80 7 

8.93 

9.39 


7.70 3 
7.96 4 




8.9 + 0.8Y 15 


0859-039 


23.78 


7.77 


0.51 


29.1 


21 


23.22 
24.20 


7.84 12 
7.88 18 






1042-690 


21.42 


7.78 


0.51 


36.3 


33 


20.64 
21.38 


7.73 6 
7.86 4 






1105-048 


15.28 


7.83 


0.52 


24.5 


142 


15.58 
15.54 
15.87 7 


7.81 4 
7.82 2 




<32 13 


2359-434 


8.66 


8.56* 


0.95 




2200 


7.76 7 
8.72 
8.85 
8.67 


8.58 4 
8.57 9 
8.83" 




«0 10 7.8±0.4Y 15 


'Bereeron etal.il 1989): : Bereeron et al. Hl992): 


; Beraeron et al. 


(2001)rBraaaaliaetal. 1 1995); ; Finlev et al. 1 1997); f Kawka et al. 


<200Ct): iKepler & NeranNl993l): tedder et al. 


1 1992); tester &Allard 


1993); 1C 


Koesteretal.N 19981): 'iKoester et aLlfeOOll): 



' TVennesN^gl) : ' IVennes et aiHl997l) . 
* gravity fixed at a value which reproduced the absolute brightness computed from the parallax. 



i.e. the brighter component is less massive. The luminosity ra- 
tio of both white dwarfs amounts to 1.5. 

WD 0346-011 (GD50): iBragaglia et al. I dl990h list this 
white dwarf as a suspected double degenerate from their ra- 
dial velocity measurements. However, the more accurate radial 
velocity measurements performed by Maxted et al. (2000) and 
from the SPY spectra indicate a constant radial velocity. GD 50 
is a very massive white dwarf (M = 1.24M ). However, our 
spectropolarimetric measurements did not detect a significant 
magnetic field. In our sample, this is the star with the broadest 
Balmer lines and, consequently, the weakest limit on the mag- 
netic field strength. 

WD 0446-789: We found a magnetic field of 4280 G in this 
relatively young white dwarf (cooling age of 21 Myr) with a 



mass (O.49M ) sli ghtly lower than the peak of the white dwarf 
mass distribution. iMaxted et al. I (2000) and our SPY observa- 
tions show that the radial velocity of this star is constant and 
thus membership in a close binary system can be ruled out. 

WD 061 2+1 77 (LTT 11818): By the time a star reaches 
~45000 K the process of H accumulation has progressed to 
the point that all white dwarfs are DA in appearance and no DB 
white dwarfs are observed in the range 30000-45000 K, the so- 
called " DB Gap". WD 0612+ 177 lies near the red edge of the 
DB Gap ( Holberg et al. 1990), where the onset of near surface 
convection can mix He into the photospheres of those stars hav- 
ing the thinnest H layers. Holbe rg et al. 1 ll 19901) reported the 
detection of a weak feature due to He i /14471 in three high 
signal-to-noise spectra of this white dwarf. From the observed 
strength of these features those authors obtained estimates of a 
homogeneous He/H ratio of log (He/H) = -2.56+0.26 or, al- 
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ternatively, a stratified H envelope mass of log (Mh/M©) ~ 
-16.6+0.3. Ho wever, subse quent optical observations of this 
WD bv iKidder etaTI d 19921) did not reveal the 4471 A fea- 
tures at the strength of the earlier observations. We also do not 
detect this feature in our two single observations of this WD. 
One possible explanation given by Kidder et al. is that photo- 
spheric He may not be distributed uniformly over the surface 
of WD 0612+177, possibly due to the presence of a magnetic 
field and a slow stellar rotation which leads to a modulation 
of the He i line strengths. Our observations suggest a 3cr up- 
per limit of the magnetic field of less than 2 kG from two dif- 
ferent observing dates. This rather tight upper limit needs to 
be taken into account when considering scenarios to produce a 
non-uniform distribution of He and H in the white dwarf atmo- 
sphere. 

WD 0839-327: The quality of our spectral fit is poor com- 
pared to other stars in our sampl e. Our result is virtu ally iden- 
tical to the value derived by Bragagli a et al. I ill 995b using the 
same technique. iBergeron et al. I J200ll) using optical and in- 
frared broadband photometry combined with a trigonometric 
parallax derived lower temperature and gravity. This might in- 
dicate that this object is an unresolved binary, consisting of two 
white dwarfs . This white dwarf is rep orted as a double degener- 
ate binary by Bragaglia et al. 1 (ll99ol) . However, the more accu- 
rate measurements done for the SPY project indicate a constant 
radial velocity ruling out a very close system. 

WD 1 042-690: This wh ite dwarf has an M type companion in 
a 0.34 days orbit llKawka et al. 120 00). Some contamination by 
line emission from the cool companion (mainly Balmer lines 
and Ca H+K) is present in our spectrum. The contaminated 
parts of the line profiles were excluded from our fits. Broad- 
band linear polarization could be measureable for this binary 
due to irradiation and reflection effects. This should, however, 
not affect our measurement of the circular polarization in spec- 
tral lines. 

WD 1 1 05-048: For this target we have discovered in the sec- 
ond of two observations the spectropolarimetric signatures 
from a magnetic field of -3900 G. WD 1105-048 resides in 
a common proper motion binary with a dM5 compa nion sepa - 
rated by 279", corresponding to 6840 AU JOswalt et ai~Fl9 88). 
The mass of this white dwarf (0.52M Q ) is very close to the peak 
of the white dwarf mass distribution. 

WD 2359-434 (L 870-2): For WD 2359-434 we have de- 
tected a magnetic field of -3100 G. Already iKoester et al. I 
(1998) reported a very shallow and narrow core of the Ha pro- 
file observed in this object. They speculated that what they see 
is only the unshifted component of a Zeeman triplet, with the 
other components shifted outside the observed spectral range 
or smeared out due to the inhomogeneity of the field. In favour 
of this speculation is the unusually high surface gravity of this 
object, since magnetic whit e dwarfs tend to be more massive 
than non-magnetic objects ( Liebertlll9 88). However, it would 



require a magnetic field > 50 kG to shift the cr components 
outside the spectral range of the Koester et al. observation. No 
further components could be detected in our spectra or the SPY 
spectra. 

We noticed during the fitting process that our FORS data 
show flat Balmer lines cores. Available UVES spectra of this 
white dwarf s how that the core o f Ha is also flat (as already 
reported bv IKoester et al. lll998l) . Since an accurate parallax 
measurement exists for WD 2359-434, we modified our fit- 
ting procedure for this object. The gravity was fixed at a value, 
which reproduced the absolute brightness of the white dwarf 
computed from the parallax. The temperature was determined 
from a fit of the line wings of the Balmer lines, excluding the 
cores. Since the absolute brightness depends on the stellar tem- 
perature as well, a few iterations were necessary, before we got 
a self-consistent solution. The derived temperature (8.6 6 kK) 
is consistent with the result of IKoeste r & Allardl dl993|). who 
fitted UV spectra taken with IUE. iKepler & Nelan | dl993h us- 
ing the same technique, derived a lower temperature (7.76 kK). 
This is possibly caused by the low value of the surface gravity 
(logg = 8.0) adopted during their fitting procedure, which is 
much lower than plausible values for WD 2359-434. 

The reason for the flat Balmer line cores of WD 2359-434 
remains a mystery. Being a white dwarf with a broad range of 
magnetic fields is a possible explanation since it would smear 
out the sigma components. However, in our polarization mea- 
surements there is no indication for such a broad range. 

6. Discussion and conclusions 

In this work we have used the spectropolarimetric capability 
of the FORS1 instrument, together with the light collecting 
power of the VLT, in order to investigate the presence of mag- 
netic fields in the range 1-10 kG for a sample of 12 white 
dwarfs. Three of the stars out of 12 normal DA white dwarfs of 
our sample, WD 0446-789, WD 1105-048 and WD 2359-434, 
exhibit magnetic fields of a few kilogauss in one or all avail- 
able observations (see Figs. 3, 4 and 5). The detection rate of 
25 % suggests now strongly that a substantial fraction of white 
dwarfs have a weak magnetic field. 

With the exception of the bright white dwarf 40 Eri B, 
for whic h a magnetic field of only 4 kG had earlier been 
detected JFabrika et al. II2003I) . WD 0446-789, WD 1 105-048 
and WD 2359-434 have the weakest magnetic fields detected 
so far in white dwarfs. In p revious exte nded searches for weak 
fields in white dwarfs (e.g. Schmidt & Smith I I 1995) only very 
few objects have been found. Only six detections have been re- 
ported for magnetic fields below lOOkG, which are not all con- 
firmed, and_OTil^_tlu^eobjects have a field weaker than about 
50 kG. lKawkaetal.1 ll2003l) have reported longitudinal mag- 
netic fields in three stars but no significant detection was made 
because their lcr error was almost as large as the observed 
value itself. They concluded that the population of white dwarfs 
with magnetic fields in excess of 1 MG is well known, but that 
lower-field white dwarfs remained undetected. 

Note, however, that our investigation is based on the av- 
eraged longitudinal component of the magnetic field, meaning 
that the maximum magnetic field at the white dwarf surface can 
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be stronger, depending on the field geometry (described e.g. by 
offset dipoles, or more complex distributions; with the underes- 
timate being larger for a more complex magnetic distribution) 
and on the orientation relative to the observer. Therefore, our 
results for the three objects with a positive detection are lower 
limits, since cancellation effects are expected. 

The population of known magneti c white dwarfs prese ntly 
comp ris es some 125 st ars ( Wickramasinghe & Ferrario 
200(j lGansickeetal.1 120021 ISchmidt et al. I 120031) . 
Wickramasinghe & Ferrario I J2000h and I Jordan I i l200ll) estab- 
lished that about 3% of all white dwarfs had a magnetic field 
above 100 kG. Results fro m the first two years of the Sloan 
Digita l Sky Survey (Gan sickeet all 120021 ISchmidt et~al~l 
2003) show the total number of known magne tic white dwarfs 
with B > 1.5 MG being 6%. Recently, luebert et al. I J2003t) 
have found that the incidence of magnetism at the level of ~ 
2 MG or greater is at least ~10%, or higher. They suggest that 
the total fraction of magnetic WDs may be substantially higher 
than 10% due to the limited spectropolarimetric analyses 
capable of detecting lower field strengths down to ~ lOkG. 
Our 3 detections out of 12 objects seem to indicate that low 
magnetic fields on white dwarfs (<10kG) are frequent while 
high magnetic fields are relatively rare. However, with only 
three detections this hypothesis remains insecure. If confirmed 
by future observations, the investigation of weak magnetic 
fields in white dwarfs could form a cornerstone for the future 
investigation of the properties and evolution of stellar magnetic 
fields. 

Our sample of white dwarfs is too small to discuss in de- 
tail the dependence of the magnetic field strength on the stellar 
parameters (masses and cooling ages). It is, however, worth- 
while to mention that two of our detections (WD 0446-789 
and WD 1105-048) have masses of only 0.5 M Q . This means 
that their progenitor s on t he main-sequence had less than 
1 M {Weidemann 2000). These two stars are therefore 
very different from the majority of white dwarfs with mega- 
gauss magnetic fields which tend to have higher masses 
jGreenstein & Oke II 1 982t iLiebert lfl98^ and, therefore, high- 
mass parent stars. 

Measurements of weak magnetic fields are now possible 
for many white dwarfs with the new large telescopes, which 
allow a magnetic field function (MFF, in analogy to the mass 
function) to be constructed in the 1-100 kG range once a suf- 
ficient number of detections have been made. Such a MFF can 
be co mpared to the corresp onding function for main-sequence 
stars JBvchkov et al. 1 19971) and will provide input for answers 
to the following key questions on the evolution of magnetic 
fields in stars: Are the magnetic fields in white dwarfs just the 
fossil relics of magnetic main-sequence stars strengthened by 
contraction due to conservation (to a large extent) of magnetic 
flux? Or do the magnetic fields develop considerably through 
the final stages of stellar evolution? Are the strongly magnetic 
white dwarfs a distinctive class of objects or do they just rep- 
resent a tail of the distribution of magnetic fields present in all 
white dwarfs? Is there a dependence between magnetic field 
strength and mass as found in the case of magnetic WDs with 
higher field strengths? Do the magnetic field strengths corre- 



late with temperature, which would be a hint for a decay on the 
white dwarf cooling sequence? 

Several authors have suggested that the frequency of mag- 
netic white dwarfs may increase with decreasing effective 



temperature, luminosity and with increasing cooling age (e.g 


IValvavin &Fabrikalll998l 


Liebert et al. 


20031), and may de 



crease sharply with distance (Fabrika & Valvavin 1998). 

Alternativ ely to the foss i l origi n of the magnetic field in 
white dwarfs. Markiel et al. I i 19941) and lThomas et all i 19951) 
have shown that a weak magnetic field of * 1.3 kG in the vari- 
able DB star GD 358 can be explained by an ato dynamo. The 
magnetic field in this star has been inferred indirectly by an- 
alyzing the g-mode oscillat ion spectrum taken with the WET 
(Whole Earth Telescope, Wing et et al. lll994T) . However, ac- 
cording to the atmospheric parameters, the convection zone in 
all of our sample stars should be too shallow to support an aco 
dynamo. 

Another long-standing problem in white dwarf research 
is the issue why metals are accreted by helium-rich white 
dwarfs in the range 8000 K < T eff < 15000 K dur- 
ing the passage through an interstellar cloud while almost 
no hydrogen is brought in to the white dwarf atmosphere. 
Illar ionov & Sunvaevl lll975) suggested that fields below 10 5 G 
provide a screening mechanism to separate hydrogen and ion- 
ized species from grains in white dwarfs accreting from the in- 
terstellar matter. Since this phenomenon always occurs in this 
type of star, it is possible tha t at this level al l white dwarfs con- 
tain magnetic fields. Friedri ch et al. I J2004f) have searched for 
circular polarisation in one DBZ and one DBAZ which have 
accreted metals, but three or four orders of magnitude less 
hydrogen than expected. In one case (L745-46A) a magnetic 
field of 7kG (1 cr error of +2kG; 99% confidence interval 
of +6kG) was found, which, however, was based on the Ha 
line only. For the second object (GD 40), only an upper limit 
of 12 kG (99% confidence) could be derived from polarization 
measurements around three spectral lines. Theoretically, mag- 
netic fields of 4 kG and 250 kG, respectively, would be required 
for the screening by the propeller mechanism to be efficient 
around the two stars. 

At our signal-to-noise ratio, magnetic fields down to about 
2 kG can be measured. There is still the possibility that all 
magnetic white dwarfs contain surface magnetic fields at the 
1 kG level. In order to test this hypothesis, much longer expo- 
sure times would be necessary, even with the VLT. 
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